Initial Stages of Bose-Einstein Condensation 
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We present the quantum theory for the nucleation of Bose-Einstein condensation in a dilute 
atomic Bose gas. This quantum theory comfirms the results of the semiclassical treatment, but has 
the important advantage that both the kinetic and coherent stages of the nucleation process can 
now be described in a unified way by a single Fokker-Planck equation. 



PACS number(s): 03.75.Fi, 67.40.-w, 32.80.Pj 

In two previous papers we have developed the semi- 
classical theory for the nucleation of Bose-Einstein con- 
densation in a weakly- interacting Bose gas At that 
time the main reason for studying this nucleation prob- 
lem, was to determine whether and on what time scale 
Bose-Einstein condensation could be observed in ongoing 
experiments with magneticaly trapped atomic hydrogen 
HI), an d f° r that purpose the semiclassical theory was 
sufficient. However, now that Bose-Einstein condensa- 
tion has actually been achieved in atomic 87 Rb Q, 23 Na 
H and 7 Li Q vapours, it appears that experimental stud- 
ies of the dynamics of Bose-Einstein condensation can be 
performed in such detail, that a more elaborate theory 
is required to fully understand the outcome of these ex- 
periments. It is the aim of this paper to present such a 
theory. 

In the semiclassical picture the formation of the con- 
densate proceeds in three stages. In the first stage the 
gas is quenched into the critical region of the phase tran- 
sition, for example by means of evaporative cooling. This 
quenching of the gas is a kinetic process and can be 
accurately described by the quantum Boltzmann equa- 
tion. However, such a kinetic equation cannot lead to a 
macroscopic occupation of the one-particle ground state. 
To achieve that, a second stage is needed in which the 
gas first develops the instability towards Bose-Einstein 
condensation and then coherently populates the ground 
state by a depletion of the low- lying excited states. Af- 
ter this coherent stage the gas has clearly a highly non- 
equilibrium energy distribution, and must come to equi- 
librium in a third and final stage. This last stage is again 
of a kinetic nature and can be studied by the appropriate 
quantum Boltzmann equation for the Bogoliubov quasi- 
particles of the Bose condensed gas. 

Although the semiclassical theory gives important in- 
sight into the physics of the nucleation of Bose-Einstein 
condensation and the various time scales involved, it is 
not completely satisfying because of the following two 
reasons. First, it makes a sharp distinction between ki- 
netic and coherent processes which in principle of course 
take place simultaneously. Second, the coherent stage of 
the evolution is described by a time-dependent nonlinear 
Schrodinger equation for the complex order parameter 



<p(x,t) = (ip(0,t)}, that formally has the property that if 
the order parameter is zero initially it will always remain 
zero. Therefore, the semiclassical theory actually makes 
use of the fact that there are quantum fluctuations in 
4>{x, t) without taking these explicitly into account. Both 
these fundamental problems are resolved in the quantum 
theory, to which we now turn our attention. 

To arrive at the quantum theory of the order parame- 
ter 4>(x, t) it is, in contrast to the semiclassical approach, 
not sufficient to derive the equation of motion for the ex- 
pectation value of the (Heisenberg) operator ijj(x, t) that 
annihilates an atom at position x and at time t. Instead 
we need to find an equation that determines the evolution 
of the full probability distribution P[<fi*, <j>; t], so that we 
are able to consider also the fluctuations in </>(x, t) . This 
is most easily achieved as follows. Introducing the initial 
density matrix p(to) of the at that time still uncondensed 
gas and the coherent states 



\cf>(x),t)=exp\ / dx4>(x)^(x,t)}\0) , 



(1) 



where |0) denotes the vacuum state, the probability dis- 
tribution P[<j>* , 4>; t] of interest equals 



; t] = Tr 



p(to) 



l0,t)(<M 
(<M<M) 



(2) 



Moreover, using an expansion of p(to) in terms of the 
above coherent states, the latter can be rewritten as the 
functional integral 



;t] = y dfoSMto 



p[\<h 



2 . iK^>*I^0j*0>I" 

>*o 



(<M<M) 



(3) 



Since (</>,i|</>,i) = exp if dx \<f>{x)\ 2 } is in fact indepen- 
dent of time, our task is essentially reduced to the deter- 
mination of the probability \(<p, t\<fio, to) | 2 . 

This we achieve by writing (</>, t\4>o, to) 'path' in- 
tegral over all field evolutions from < to t. Similarly, 
we also write (<p,t\<po,to)* — (4>o,t \(j>,t) as a 'path' in- 
tegral over fields, but these can now be seen as evolving 
backwards in time from t to to- Performing then also the 
integrations over 4>o (x) in Eq. (^) , we see that we are lead 
in a natural way to the Keldysh formalism developed in 
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Ref. [pj, because we have to perform a functional inte- 
gral over all field configurations ip(S,t) that evolve from 
to to t and back to to, i-e. along the Keldysh contour C. 
Using the results obtained there, we conclude that the 
probability distribution for 4>{x) is given by 



p[<t>*,4>;t] = i dm<m ex P <j -s\r^] \ , (4) 



where the effective action S[ip*,ip] in the integrand be- 
comes 



SWM = E {^(*) + ^(*Ws - l^l 2 } + j£ dtj dt' 

k 

~w E ^(*^.^;*.0^ / ^(*)^ /2 _ jS (*)^ /2 _ & (* , )^ /a+iE ,(*') 



k,k',K 



(5) 



introducing the notation = h 2 k 2 /2m — /xq = — /zo for the kinetic energy of an atom with momentum Tifc 
and mass m relative to the chemical potential /zo at ioj &{t,t') for the (5- function on the Keldysh contour defined by 
J c dt 1 S(t, t') — 1, and U for the volume of the system. 

We arrive at this result by making use of the so-called ladder or many-body T-matrix approximation. Physically, 
this means that we include all two-body scattering processes in the gas but neglect three (or more) particle processes. 
This is justified for a weakly- interacting Bose gas because of the smallness of the gas parameter na 3 , where n is the 
density and a the two-body scattering length of the atoms in the gas. Furthermore, in the same approximation the 
selfenergy obeys the Hartree-Fock like expression 



hE(k;t,t') = 2i 



dk' 



V(k -k',k- k', k + k'; t, t')G{k'; t', t) , 



(6) 



and the effective interaction between the atoms is V(k, k', K; t, t') = (T(k, k', K; t, t') + T(—k, k' , K; t, t'))/2 in terms 
of the many-body T-matrix that is determined by the Lippman-Schwinger equation 



T(fc, k', K- 1, t') = V(k - k')S(t, t') 

i r r dk» 



c 



(2tt) 3 



V(k - k")G(K/2 + %"■ t, t")G{K/2 - k"; t, t")T(k", k' , K; t" , t') 



(7) 



with V(k — k') the Fourier transform of the interatomic 
potential, 

G(x,t;x',t') = -i Tr [p(t )T c (^(x,t)^ (x' ,t'))] (8) 

the one-particle Green's function and Tc the time- 
ordening operator along the Keldysh contour. 

To extract the equation of motion for P[(j>* ,(j>;t] from 
Eq. (Q), we want to express the functional integral over 
the fields ip(x,t + ) and ip(x,t-), which respectively live 
on the forward and backward branch of the Keldysh con- 
tour, as a functional integral over fields that are de- 
fined on the real axis. This is convenient, because in 
this manner we can from the result immediately read 
off the appropriate 'Schrodinger' equation for the prob- 
ability distribution. We achieve this goal in the follow- 
ing way: First, we perform the variable transformation 
ip(x,t±) = cj)(x,t) ± x(x,t)/2 in Eq. (§]). Substituting 
this transformation into the effective action S[tp*,ip], it 
turns out to be convenient to decompose the various 
functions F(t,t') on the Keldysh contour into their re- 
tarded, advanced and Keldysh components (t, t') 
and F K (t,t'), respectively 0. Second, we notice that 
for the ultracold atomic gases of interest the thermal de 



Broglie wavelength A = {2ttT% 2 /mkBT) 1 / 2 is so large com- 
pared to the scattering length a that we are allowed to 
make a Markovian approximation to the many-body T- 
matix and hence to the selfenergy hT,(k; t, t'), i.e. we can 
use F^^ K (t,t') ~ F^' K (t)S(t-t^) for these quantities. 
Third, we expand the action up to quadratic order in 
the field xi^^), which physically describes the quantum 
fluctuations in the order parameter (ip(x,t)} that are ne- 
glected in the semiclassical theory, and then integrate out 
this field. It is important to note that the cubic terms 
in the action should indeed be neglected here to avoid a 
double counting of the effects of the interaction Q . 

After performing the integration over x(x, t) we clearly 
arrive at an expression for P[<f>*, 4>; t] as a functional in- 
tegral over (f>(x,t). Using for example the methods of 
canonical quantization, we are now in a position to derive 
the 'Schrodinger' equation that is solved by this func- 
tional integral. Carrying out the above procedure re- 
quires some straightforward but tedious algebra and ul- 
timately leads to the central result of the present paper, 
i.e. a Fokker-Planck equation describing both the kinetic 
and the coherent stages of the condensation process. It 
reads M 
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k k \ k',k" 
k & \ 

-5 E t^t (W + 1 E ^ & % 



'k'+k"-k^k"^k' 



p 



p 



(9) 



fc.fc' K 

where V (±) ' K {k, k', k") is a shorthand notation for 
V&> K {k- (k' + k")/2, (k' - k")/2, k' + k"; e(fc') + e(fc")) 
and implicitly also depends on P[4>* ,4>]t], because it 
depends on the average occupation numbers N{k\ t) = 
Ti[p(to)ipt(t)ipj:(t)}. The same is in fact true for the var- 
ious selfenergies as we will see now in more detail. 

In the above Fokker-Planck equation the nonlinear 
terms proportional to V^' K (k, k', k") should be re- 
garded as giving corrections to the linear terms propor- 
tional to the selfenergies KS^ ,K (k; t) due to the con- 
densate. The effects of the above condensate particles 
have already been included correctly in the linear terms 
and should not be accounted for twice. Therefore, in 
the normal state of the gas we must use the linearized 



version of Eq. (^|). Restricting ourself to this case first, 
we see that the real part of HJ^^ (k;t), which we call 
S(k; t) in the following, describes the change in the in- 
stantaneous energy levels of the gas due to the inter- 
actions. Indeed, the instantaneous energy levels obey 
Tiw{k]t) = e(k) + S(k;t). From our experience with the 
semiclassical theory we know that the selfenergies change 
on a very short time scale of 0(h/kBT c ) for conditions 
near the critical temperature T c . However, the change of 
the occupation numbers, as derived from the imaginary 
part R(k;t) of h'S^ (k-.t), occurs only on a much longer 
time scale of 0((A c /a) 2 Ti/kBT c ) determined by the mean 
free path of the atoms. As a result it is very accurate to 
apply an adiabatic approximation and use 



S(k;t) ~ 2nT (+) (0,0;0) 



dk' 



dk" 



V {+) {0,k',k") N(k';t)N(k";t) 



1 - cos(7jfc' • k"(t - t )/m) 



(2tt)3 J (2tt)3 ' ' ----- -v- tfk'.k'i/m 

where TW(0,0; 0) = 4nah 2 / m is the two-body T-matrix at zero momenta and energy, and 
dk' f dk" 



(10) 



V {+) (k, k', k") 5 ( e(fc') + e(k") - e(k' + k" - k) - e{k) 



Rik ^ = -^J {27T) 3 J (27r) 3 

(1 + N(k'; t))(l + N(k"; t))N(k' + k" -k;t)- N(k'; t)N(k"; t)(l + N(k' + k" - k; t)) 

as also expected from a Fermi's Golden Rule calculation of dN(k;t)/dt. Moreover, in the same aproximation 



(11) 



dk' 



dk" 



V W (fc, k', k")\ 2 5 (e(fc') + e{k") - e(k' + k" - jfe) - e(jfe)) 



(2tt) 3 J (2tt) 3 

(1 + N(k'; t))(l + N(k"; t))N(k' + k" -k;t) + N(k'; t)N(k"; i)(l + N(k' + k" - k; t)) 



(12) 



Notice that in equilibrium the above relations imply that 
h^ K (k) = 2i(l + 2N(k))R(k), which is just the famous 
fluctuation-dissipation theorem. 

In the normal state the Fokker-Planck equation in 
Eq. (^) thus combines the kinetic evolution of the occu- 
pation numbers due to the quantum Boltzmann equation 
with the coherent evolution of the energy levels given by 
Tiuj(k;t) ~ e(k) + 5(0; t). In particular, the solution of 
the Fokker-Planck equation in the limit t — > oo is due to 
the fluctuation-dissipation theorem correctly given by 

P[0*,0;oo] = ]Ja(fc)exp|-a(fc) |</) E | 2 | , 



with a(k) = l/(N(k) + 1/2), N(k) the Bose distribu- 
tion evaluated at e(fc) + 5(0; oo) — /j,q and 5(0; oo) deter- 
mined selfconsistently by Eq. (|l0|). In equilibrium, Bose- 
Einstein condensation therefore occurs at a temperature 
such that 5(0; oo) = /io. From the semiclassical analysis 
we know that this can only be fulfilled in the case that 
a > and then leads to a critical temperature T c for 
the interacting gas that is higher than that for the non- 
interacting gas. Interestingly, the same conclusion was 
recently also reached by a completely different calcula- 
tion || . 

Returning to the nonequilibrium problem of inter- 
est, we see that if the gas is quenched sufficiently far 
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into the critical region, the effective chemical potential 
(icl(t) = A*o — 3(0; t) will become positive during the ki- 
netic evolution of the gas towards equilibrium. At that 
time, say at t = t%, the gas is unstable against the forma- 
tion of a condensate. To determine how the condensate 
actually grows after the instability has occured, we must 
return to our full nonlinear Fokker-Planck equation and 
discuss the coherent stage that now follows. Neglecting 
therefore the kinetic, i.e. imaginary, terms in the right- 
hand side of Eq. (^|) and making use of the fact that in 
the coherent stage of the evolution we are mainly inter- 
ested in momenta hk < hy/na, we can easily show that 
a solution to the Fokker-Planck equation is 



;t] = J d^lWi 



pm 2 -Mm 2 -\<t>d\ 



where P[\<t>i \ 2 ',ti] is the probability distribution after the 
kinetic stage and cf> c i (x, t) is the solution of the nonlinear 
Schrodinger equation 



th- 



at 



2m 



/i c; (i)+T(+>(0,0,0;0)|</> c; | 2 U d 



with the initial condition (f) c i(x,t\) = <j>\{x). 

Clearly, with this approximate solution to the Fokker- 
Planck equation we have essentially rederived the semi- 
classical theory of the coherent stage. It is interesting 
to note that in contrast with the quantum theory for 
the turning on of a laser |l(J , this coherent stage of the 
evolution is important. The reason for this is that in 
the case of a laser there is only one mode for the pho- 
ton field and the nonlinear Schrodinger equation then 
only affects the phase of this mode, with the result that 
P[(j)* ,4>;t] — P[\(j>\ 2 ;ti\. However, in a Bose gas we are 
dealing with many modes and the nonlinear Schrodinger 
equation actually leads to the building up of population 
in the one-particle ground state due to a depopulation of 
the lowest excited states. Introducing the phase \ci of 
/ dx <j) c i, the latter picture arises roughly speaking be- 
cause the instantaneous energies of the excited states are 
hw(k;t) = [(e(fc) + ?ixd-M c / + 2n c T(+)) 2 -KT(+)) 2 ] 1 /2 j 

which are imaginary for low momenta since the time 
derivative h Xc i(t) < fi d (t) - n c (t)T(+'(0, 0, 0; 0) as the 
population of the lowest excited states decreases in favor 
of the condensate density n c (t). 

Having formed a small condensate density in this way, 
the gas then has to evolve towards an equilibrium be- 
tween the condensate and the noncondensed part of the 
gas. For this the imaginary or kinetic terms in the right- 
hand side of Eq. (^) are responsible. Indeed, using the 
explicit expressions for V^' K (k,k' ,k") following from 
Eq. (Q), we find that in lowest order the presence of the 
condensate has the effect of replacing N(k' + k" — k; t) 
by N(k' + k" -k;t) + Vn c (t)(2Tr) 3 5(k' + k" - k) in the 
normal state expressions for R(k;t) and hY; K (k;t). The 



kinetic equations for the superfluid state of the gas are 
therefore essentially the same as those recently studied 
by Semikoz and Tkachev Although these do not 

take properly account of the coherence effects at low mo- 
menta hk < hy/na, they are sufficiently accurate for our 
purposes because for the quantum gases of interest we 
always have that fi/A> hy/na near the critical temper- 
ature. 

This completes our discussion of the quantum theory 
for the formation of a condensate in a weakly-interacting 
Bose gas. Summarizing, we have derived a Fokker-Planck 
equation for the probability distribution of the order pa- 
rameter (ip(x, i)), which gives a desciption of both the 
kinetic and coherent stages of Bose-Einstein condensa- 
tion. This Fokker-Planck equation implies that the order 
parameter obeys a nonlinear Schrodinger equation with a 
Langevin noise term tj(x, t), which in a good approxima- 
tion has a Gaussian probability distribution that fulfills 
the fluctuation-dissipation theorem. We have also dis- 
cussed the qualitative behavior of the solutions to the 
Fokker-Planck equation. We hope to report on a more 
quantitative numerical study in the near future. 
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